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^ ' Abstract 

o: 

^' . We present a class of exactly solvable quantum spin models which consist of two 

Heisenberg-subsystems coupled via a long-range Lieb-Mattis interaction. The total 

system is exactly solvable whenever the individual subsystems are solvable and allows 

to study the effects of frustration. We consider (i) the antiferromagnetic linear chain 

^ . and (ii) the Lieb-Mattis antiferromagnet for the subsystem-Hamiltonians and present 

10 ' (i) the complete ground-state phase diagram and (ii) the full thermodynamic phase 

f^ . diagram. We find a novel phase which exhibits order from disorder phenomena. 



PACS: 75.10, 75.10J, 75.50E 



Introduction. ~ The properties of interacting quantum spins have been studied inten- 
sively over a long period. Some selected model Hamiltonians, which typically describe 
homogeneous collections of interacting spins, can be solved exactly. Examples are (i) 
models in one dimension solvable by Bethe-Ansatz [Q , (ii) a one-dimensional model with 
long-range inverse-square exchange [0, (iii) valence-bond models 0] and (iv) models with 
long-range interaction of constant magnitude (Lieb-Mattis type models 0| which have 
recently been used to discuss spontaneous symmetry breaking in spin systems Q). Our 
understanding of interacting quantum spins has benefited greatly from these studies. 

On the other hand, in experimental applications we have to deal sometimes with 
situations where two subsystems of spins interact mutually. One example is the garnet, 
Mn3Cr2Ge30i2, in which the Mn and Cr moments are located on two disjunct, inter- 
penetrating sublattices |y, 0, with the intra-sublattice couplings dominating over the 
inter-sublattice coupling due to frustration effects [^, |8|. Another example is the two- 
layer structure in certain high-Tc compounds, like YBa2Cu306 P], where the competition 



between the inter-layer and the intra-layer couplings might lead to a spin-gap [10|. 

Here we present a class of exactly solvable models which consist of two subsystems 
coupled via a long-range Lieb-Mattis Q interaction. The total system is exactly solvable 
whenever the solutions of the individual subsystems are known and it allows to investigate 
the effects of competing interactions. We consider two different types of subsystems, 
the linear chain and the Lieb-Mattis model, and present the ground-state (GS) phase 
diagrams. In addition we present for the case of two coupled Lieb-Mattis subsystems 
the full thermodynamic phase diagram. We find a novel phase which exhibits order from 
disorder phenomena fl^, S] . 

Model. - We consider a system of A^ = A^i + A2 interacting quantum spins {s = 1/2) 
made up of two mutually interacting subsets described by the Hamiltonian 

H = aiHi{si} + a2H2{si} + aHint{si,si} (1) 

with a, ai,a2 > 0. Here Hi{si} (i = 1, . . . , Ni) and H2{si} (/ = Ai + 1, . . . , A^) contain 
the interactions within the first and the second subset of Ai and A2 spins respectively. 
The antiferromagnetic (a > 0) interaction between the two subsets of spins is taken to 
be of the Lieb-Mattis form ||, Hint{si,si} = Si • S2/(AiA2) with Si = Eil\ Sj and 
S2 = J2i=Ni+i ^^ We assume Hi and H2 to be spin-rotationally invariant. Then the oper- 
ators S^ = (Si + S2)^, Sf , S2, H , Hi and H2 all commute with each other. The eigenval- 
ues of Eq. (|) are given by Es,s^,s, = [a/{2NiN2)] [S{S + 1) - Si{Si + 1) - 52(^2 + 1)] + 
aiEi + Q2-E2; where Si, S2, Ei, E2 and S denote the quantum numbers of Sf, S2, Hi, H2 
and S^ respectively with S £ [|5'2 — ^il, ^2 + Si]. We consider the thermodynamic limit, 
Ai, A2 — > 00, and define the normalized quantum numbers x-y = 2S^/N^ (0 < x-y < 1, 
7 = 1,2). The GS energy which is realized for S = 152 — ^il, can then be written as 
Eq = —a X1X2/4 + ai-Ei(xi) + a2E2{x2), where Ei{xi) and E2{x2) are the lowest eigen- 
values of Hi and H2 within the subspace of a given magnetisation xi and X2. For a large 
class of antiferromagnetic Hamiltonians H^ a level ordering E^{x) > E^(x') is valid for 
x > x' [^. Then a competition arises in Eq between the term proportional to a which 
tends to maximize x^ and the terms proportional to a^ which tend to minimize x^, giving 
rise to frustration effects. 

We consider two models for the Heisenberg Hamiltonians Hi{si} and H2{si} entering 
Eq. (H) and set a = 1 for the rest of this paper. The two models are the antiferro- 
magnetic linear chain (LC), HIj''-^ = (1/A^) ^^ SjSj+i, and the Lieb-Mattis (LM) model, 
if^^ = (2/A'^)2S^S;^, with S:^ = E^eAe'y^i and S^ = EieBey^i being the total-spin 



operators of the A, B sublattice of the respective subsystem (7 = 1,2). The scaHngs 
with A'^ are chosen such that (i) ah energies become intensive in the thermodynamic 
Umit A'^i, A'^2 -^ 00 and (ii) that for the fuhy polarized state {x-y = 1) the energies are 
£'^C'(1) = E^'^''^{1) = 1/4. Both models can be solved exactly. For the LM model in 
addition to S^ also the respective sublattice spin-operators, (S^'-^)^, commute with the 
Hamiltonian and the corresponding (normalized) quantum numbers, x^'^ = {4/N^)S;^' 
(0 < x^'^ < 1) are separately conserved. 

For the LM model, mean-field theory becomes exact and quantum fluctuations van- 
ish identically in the thermodynamic limit. The energy E^{x^) can be calculated as 
Ei^ = 2;^/2 — -g. On the other side, quantum fluctuations are maximal for the antiferro- 
magnetic LC. We computed Ei" (x^) by numerically solving the Bethe-Ansatz equations 
for finite magnetisation Xj Q. 

Ground-state phase diagram. - The knowledge of £'^(x^) (7 = 1, 2) allows to determine 
the GS phase diagram, which we present in Fig.l. The phase diagram is symmetric 
under ai <-> 02 and the GS has minimal total spin S = {82 — Si\. For small ai and 
02 the antiferromagnetic inter-subsystem interaction dominates and the subsystems are 
saturated ferromagnets (xi = X2 = 1) with parallel sublattice magnetisations. For large ai 
and 02 the antiferromagnetic intra-subsystem interaction dominates and two subsystems 
effectively decouple into their individual singlets {xi = X2 = 0)- A line of first-order phase 
transitions separates the spin-singlet region from a region near the axis where one of the 
subsystems is partially magnetized and the other one is fully magnetized. The states with 
partial magnetisation are connected to the fully polarized state by a second-order phase 
transition. 

In Fig. 2 we have plotted the GS expectation values of the inter-subsystem correlation 
function, (sjS;)i^2 (i and I are in the first and second subsystem, respectively), and the 
intra-subsystem correlation function (sjSj)i^i of the first subsystem (for the LC case j = 
i + 1, for the LM case i £ A and j £ B) as a function of ai for fixed 02 = 0.075. For small 
ai saturated ferromagnetic long-range order within the subsystems ((sjSj)i^i = +1/4 but 
(sjS;)i^2 = —1/4) is present up to af = 1/4 for both the LC and the LM case. At this 
point a second-order phase transition leads to the non-saturated magnetic structure. For 
the case of two coupled LM systems we have that (sjSj)i^i — 1/4 ~ (ai — 1/4) while for 
the case of two coupled LC's (siSj)i^i — 1/4 ~ ^/ai — 1/4, both for small (ai — 1/4). The 
latter result can be rigorously established from the analytic expression Ei"'-^{xy) — 1/4 = 
— (1 — Xy) + (7r^/48)(l — x-yY + 0((1 — x^)^) for the energy of the antiferromagnetic LC 
near maximal polarization, x-y, which can be extracted from the Bethe-Ansatz |l|. 

The first- and second-order phase boundaries shown in Fig.l meet in tricritical points 
at (01,02) = 1/4 (1 , — 1 + jj^) and at (01,02) = 1/4 (—1 + j^ , 1) for the LC case and 
in a tetracritical point at 1/4 (1, 1) for the LM case. The line of the first order transition 
is the result of a numerical solution of the Bethe-Ansatz for the LC system, whereas for 
two coupled LM systems it is given by 160102 = 1- Finally we note that the straight first- 
order line between the two tricritical points in Fig. 1(a) separates the state with classical 
ferromagnetic long-range order within the subsystems (x-y = 1,7 = 1,2) from the quantum 
disordered state without long-range order (Bethe singlet). 

Thermodynamics. - For the Lieb-Mattis model the saddle-point approximation be- 
comes exact, i.e. the partition function is given in the thermodynamic limit by its largest 
term [11|. The same statement holds also for the model of coupled Lieb-Mattis subsys- 



tems considered here and all thermodynamic potentials can be determined. The quantum 



numbers of the largest term in the partition function depend on both ai, 02 and the tem- 
perature T. For the Lieb-Mattis model the four (normalized) quantum numbers x^'^ of 
the total spin operators (S^'^)^ of the respective sublattices are independently conserved. 
Due to symmetry we have x^ = x^ = y^ (7 = 1,2) and the total subsystem magnetisa- 
tion x-y G [0, y-y]. The subsystem magnetisations x^ = 2S^/N^ are always antiparallel and 
5 = 15*2 — Sil- We can define three independent order parameters V'l, '02 and ^3 by 

V'l = 2/1 - a;? ; V2 = y^ - ^^l ; il^s = xf + xj . (2) 

The parameters Vi and ip2 are proportional to the antiferromagnetic order parameter, 
ip^ = 4{1/N^){{S^ — S^y), of the respective subsystems (7 = 1, 2) and ips is proportional 
to the order parameter for ferromagnetism in the individual subsystems, 
V'a = 4[(l/7Vf)(Sf) + {l/N^){Sl)] . Though there are 2^ = 8 possible phases, one of 
them, ^j > 0, i = 1, 2, 3, is found to be thermodynamically never stable. In the table we 
list the stable phases and illustrate in Fig. 3 the phase diagram as a function of ai and 
temperature T for a constant 0:2 = 0.075. The high-temperature paramagnetic P-phase 
has the highest symmetry. For small a^ (7 = 1,2) the inter-sublattice exchange domi- 
nates and both the first and the second subsystem are ferromagnetic. We call the resulting 
phase the ferromagnetic, or F-phase. For large a^ (7 = 1, 2) the intra-sublattice couplings 
dominate, leading to a singlet state in both subsystems, corresponding to the collinear 
antiferromagnetic AF-phase. The low temperature AF-phase gives way to the AFi-phase 
at higher temperature when the second subsystem becomes paramagnetic, X2 = (the 
AF2 phase is realized for certain ai < 02)- 

The symmetry group of the noncollinear "order-from-disorder" ODi-phase (the OD2- 
phase is realized for certain ai < 02) is a subgroup of the symmetry groups of both 
the F-phase and AFi-phase (compare the table), which in turn are subgroups of the 
fully symmetric P-phase. These relations are consistent with general requirements for the 



tetracritical point given by ai = J{l/4)a2 + 1/8 — (l/2)a2 and ksT = (l/2)ai, in which 



these four phases meet [12|. All phase transitions shown in Fig. 3 are of second order (with 
the characteristic jump in the specific heat), besides the transition from the ODi-phase 
to the AF-phase. The symmetry group of neither the ODi-phase nor of the AF-phase is 
a subgroup of the other (see the table) and the transition is consequently of first-order. 

For every region of the phase diagram a distinct set of self-consistency equations deter- 
mines the free energy and the correlation functions. E.g., the subsystem magnetisations 
xi = 2:2/(401) and X2 = y2 as well as the sublattice magnetisations of the first subsystem, 
yi, in this ODi-phase are determined via the equations 

yi = ta,VLh.[aiyi/{2kBT)\ ; 2:2 = tanh [ (l/oi — 802) X2/(16A;b?') ] • (3) 

The inter- and intra-subsystem correlation functions for the ODi-phase are given by 

(siS;)i,2=-x|/(16ai) ; (s,s„)2,2 = x^ /4 ; (s,Sj-)i,i=2;| / (32af ) - yf /4 , (4) 

where i,l £ A and j,m £ B in (sjSj)i^i and (8/8^)2,2- The transition from the ODi-phase 
to the AFi-phase takes place when X2 = in Eq. (^) and the transition to the F-phase 
takes place when yi = xi = X2/(4qi). 

In the ODi-phase the sublattice magnetisations in the first subsystem are neither par- 
allel nor anti-parallel, one might call the ODi-phase also a "twisf'-phase. The temperature 
dependence of the inter-sublattice correlation functions in the first subsystem, (sjSj)i^i, 



are presented in Fig.4 for some selected values of ai and 02 = 0.075. In the ODi-phase 
(1/4 < ai < 1/(1602) at r = 0) the correlation function (siSj)i^i increases in magnitude 
as a function of temperature. This effect is particulary pronounced for ai = l/v8- For 
ai = 1/ v8 the A- and the B- sublattices of the first subsystem are completely uncorrelated 
at zero temperature (see Eq. (Q)) and become more and more correlated with increasing 
temperature, as long as we remain in the ODi-phase. This phenomenon is called order 



from disorder |13,|g], and is due to the interplay between different competing energy scales 
govering the correlations (compare Eq. (y) and Eq. (^)) of the individual subsystems. The 
disappearance of {siSj)i^i at zero temperature for ai = l/yo and the increase of magnetic 
order by fluctuations is related to the maximum in the GS energy Eo{ai) which signals 
maximal frustration. On the other hand, the correlations decrease in the usual way with 
increasing temperature within the ferromagnetic phase or the antiferromagnetic phase. 

Summary. - We have discussed quantum spin Heisenberg systems built of two coupled 
subsystems. If the individual subsystems are integrable on their part, the whole system 
is solvable for a long-ranged inter-subsystem coupling of constant magnitude. The three 
antiferromagnetic exchange integrals, namely the two intra-subsystem couplings and the 
inter-subsystem exchange parameter compete with each other. This competition gives 
rise to interesting frustration effects. Varying the strength of the three couplings we can 
tune the magnetic ordering of the individual subsystems between ferromagnetic and anti- 
ferromagnetic. In the limit of maximal competition between different exchange couplings 
neither ferromagnetic nor antiferromagnetic correlations dominate in the ground state. 
However, additional thermal fluctuations change the interplay between the competing en- 
ergy scales and lead to an increase of magnetic order with temperature. This mechanism 
driving the order from disorder scenario seems to be not restricted on the special Lieb- 
Mattis coupling and we argue that the same phenomenon should be present in systems 
with more realistic couplings. 

We are indebted to N.B. Ivanov for fruitful discussions. This work was supported by 
the DFG (Ri 615/1-1). 
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Table 1: The seven thermodynamically stable phases for two coupled Lieb-Mattis antifer- 
romagnets. y^ = {4/Nj)S^'^ (0 < y-y < 1, 7 = 1,2) are the relative magnetisations of the 
A,B-sublattice of the first/second subsystem and x^ = {2/N^)S^ (0 < y-y < 1) x^, G [0,y-y] 
are the magnetisations of the respective (7 = 1,2) subsystem. The paramagnetic (P), fer- 
romagnetic (F), the three antiferromagnetic phases (AFi), (AF2), (AF) and the order from 
disorder phases (ODi) and (OD2) are characterized by the order parameters t/ji = yl 
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Figure 1: The ground-state phase diagram of a) two coupled antiferromagnetic linear 
chains and b) two coupled Lieb-Mattis antiferromagnets, as a function of the intra- 
subsystem couplings ai and 02- The inter-subsystem coupling is a = 1. The phases 
are characterized by the magnetisations xi and X2 of the subsystems. The dashed/full 
lines denote phase transitions of first/second order. 



Figure 2: The ground-state inter-subsystem (sjSi)i^2 and the intra-subsystem (sjSj)i^i 
correlation function for two coupled linear chains (LC, dashed lines) and two coupled 
Lieb-Mattis antiferromagnets (LM, full lines), as a function of ai. The intra-subsystem 
coupling a2 = 0.075 and the inter-subsystem coupling a = 1. At ai = 1/4 a phase 
transition of second-order occurs for both systems and the correlation functions have a 
square-root singularity (kink) for the case of two LC's (LM's). 



Figure 3: Phase diagram of two coupled Lieb-Mattis antiferromagnets, as a function 
of temperature T and the intra-subsystem coupling ai, for 02 = 0.075 and the inter- 
subsystem coupling a = 1. The dashed/full lines denote phase transitions of first/second 
order. The phases are labelled with respect to the values of the correlation functions 
within the first subsystem. 



Figure 4: The intra-subsystem correlation function (sjSj)i^i (z G ^, j G B) for two coupled 
Lieb-Mattis antiferromagnets, as a function of temperature T for some selected ai's, and 
fixed 02 = 0.075, a = 1. The kinks in (sjSj)i^i correspond to phase transitions of second 
order. Note the particular strong "order from disorder" phenomenon for ai = 8~^'^. 



